Theory of the asymmetric ripple phase in achiral lipid membranes 
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We present a phenomenological theory of phase transitions in achiral lipid membranes in terms of 
two coupled order parameters - a scalar order parameter describing lipid chain melting, and a vector 
order parameter describing the tilt of the hydrocarbon chains below the chain-melting transition. 
Existing theoretical models fail to account for all the observed features of the phase diagram, in 
particular the detailed microstructure of the asymmetric ripple phase lying between the fluid and 
the tilted gel phase. In contrast, our two-component theory reproduces all the salient structural 
features of the ripple phase, providing a unified description of the phase diagram and microstructure. 

PACS numbers: 87.16.D-,61.30.Dk 



Phospholipids self-assemble in water to form a rich 
variety of spatially modulated phases [l| . The simplest 
of these is the 1-dimensionally modulated fluid lamellar 
phase (L a ) consisting of periodic stacks of lipid bilayer 
membranes separated by water, where the hydrocar- 
bon chains are floppy with liquid-like in-plane order. 
Changing the temperature or water content induces a 
sequence of symmetry breaking transitions character- 
ized by unique microstructures. 

On reducing the temperature below the chain melting 
(main) transition (T m ), the L a phase of phosphatidyl- 
cholines (PCs) transforms to a gel phase (Lpi), charac- 
terized by fully-stretched all — trans chains which are 
tilted with respect to the bilayer normal In ad- 

dition, an asymmetric ripple phase (Pp>) is found to 
occur in between the L a and L g > phases in many PCs 



at high water content 



Ilia. 



Extensive studies using a variety of experimental 
techniques 0, [B"18|, reveal that the P a i phase is char- 
acterized by a periodic saw-tooth height modulation 
of the bilayers having an amplitude of ~ 1 nm and a 
wavelength of ~ 15 nm, and a bilayer thickness that is 
different in the two arms of the ripple (fig. 
As a result, the rippled bilayers lack a mirror plane nor- 
mal to the rippling direction. While in principle, this 
discrete symmetry breaking can arise from an asymme- 
try in either shape (unequal lengths of the two arms) or 
bilayer thickness (unequal bilayer thickness in the two 
arms), in practice these asymmetries seem to appear 
simultaneously. 

At first it was believed that the origins of the asym- 
metric ripple lay in the chirality of lipid molecules [20j . 
However, subsequent experiments usin g ra cemic mix- 
tures showed this was not the case Il5j. More re- 
cently, all-atom molecular dynamics simulations of lipid 
bilayers have observed that the degree of chain order- 
ing is different in the two arms of the ripple |2l| . The 



occurrence of the ripple phase only in those lipids that 
exhibit a Lai phase at lower temperatures [22j, and in 
isolated bilayers |23|, suggests an intimate connection 
between chain tilt and the ability of the bilayers to form 
ripples. 



Several theoretical models have been proposed to de- 
scribe the sequence of phase transitions in such lipid 
bilayers and the microstructure of the ripple phase 
[20I . 24 - 34 1 . None of them accounts for all the ob- 
servations. We list three key features that should be 
explained by any theory of the ripple phase in achiral 
bilayers : (1) occurrence of phase between L a and 
Lpi phases, separated by two first-order transitions; (2) 
unequal bilayer thickness in the two arms of the ripple; 
and (3) unequal lengths of the two arms. 

In this paper we present a phenomenological Landau 
theory to describe the ripple phase in an isolated, achi- 
ral lipid bilayer. Our free energy expression is written 
in terms of two order parameters: a scalar order pa- 
rameter if> and a 2-D vector order parameter m (fig[2]). 
if describes the melting of the bilayer and is the dif- 
ference in the bilayer thickness [28|, [29[ between the 
fluid (L a ) and ordered (Pg/ and Lgi) phases. Since the 




FIG. 1: Electron density map of the ripple phase of dimyris- 
toylphosphatidylcholine calculated from x-ray diffraction 
data [19J. Bands labeled h and c correspond to the head- 
group and hydrocarbon chain regions of the bilayer; w de- 
notes the water layer separating the bilayers. 
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bilayer thickness is determined by the conformations 
of the hydrocarbon chains, ip can also be interpreted 
in terms of differences in chain conformations between 
the fluid and ordered phases [30(. m is the projection 
of the molecular axis on the bilayer plane |2Cj . A third 
order parameter h, describing the height of the bilayer, 
can be integrated out of the expression for the total free 
energy density. This model is found to capture all the 
three salient features of the ripple phase listed above. 



The tilt free energy density can be written as, 




FIG. 2: The unit vector n represents the orientation of the 
long axis of the lipid molecules relative to the bilayer normal 
N. m = n (N ■ n) N is the projection of n on the bilayer 
plane. 

The total free energy per unit area is taken to be the 
sum of three terms; the stretching free energy density 
f s , the tilt free energy density ft, and the curvature 
free energy density f c . For an isolated lipid bilayer f s 
is given by [Hj], 
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where if>(x,y) — S ^ x ' V g ~ S " , S(x,y) being the mem- 
brane thickness at position (x, y) in the bilayer mea- 
sured with respect to a flat reference plane, and So the 
constant thickness of the membrane in the L a phase. 
if) is taken to be positive for T < T m due to the stretch- 
ing of the chains. This is valid in general, even if the 
chains are tilted below T m . Explicit temperature de- 
pendence is assumed to reside solely in the coefficient 
of if) 2 : a 2 — a' 2 (T — T*), T* being a reference tempera- 
ture. 013 is taken to be negative, so that the continuous 
transition at T* is preempted by a first order melting 
transition at T m = T* + ^r^L. The coefficient C 
can cither be positive or negative, but 0,4, D and E are 
always positive to ensure stability. With C > 0, the 
equilibrium phases are always homogeneous in space; 
either as L a or Lp (Lai). However, with C < modu- 
lated phases are possible with some characteristic wave 
vector qo- The (V'0) 4 term is included, since in the 
context of a one dimensional model with a scalar order 
parameter, it has been shown that such a term is nec- 
essary to stabilize a modulated phase with a non zero 
mean value of the order parameter 35| . 
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i) 2 + r 3 (V-m) 4 + r 4 ^|m| 2 



+ r 5 (m • wo 2 + r 6 (m x w) 2 

+ r 7 (vV) (V • m) 2 + r 8 (v^) 2 (V • m) 2 (2) 

The first four terms in eqn.([2]) are the usual terms in 
the expansion of the free energy in terms of a vector 
order parameter. The (V ■ m) 4 term is included to be 
consistent with the (V"0) term introduced in f s . The 
next term represents the coupling between ijj and m, 
which is responsible for the appearance of tilted phases 
in this model, as 62 is taken to be positive. If > 0, 
the stable phase below T m is Lp with |m|=0. On the 
other hand, tilted phases can form if T 4 < 0. The suc- 
ceeding two terms take into account the anisotropy of 
the tilted bilayer. The next two terms represent higher 
order couplings between modulations in ip and in m, 
allowed by the symmetry of the system. 




FIG. 3: Phase diagram in the T — C plane calculated from 
the model. a' 2 = 159.42 k B , T* = 260.0 K. Values of the 
other coefficients in units of ksT* are: 03 = -306.5, 04 = 
613.15, b 2 = 0.2, 64 = 200.0, D = 557.41, E = 600.0, Tj = 

0.010, r 2 = 1.8O, r 3 = 500.0, r 4 = -3.0, r 5 = -20.0, r 6 = 

-20.0, F 7 = -500.0, r 8 = -750.0. Both the main-transition 



first order 



Ppi) and pre-transition {Ppi 



are 



The curvature energy density of the bilayer can be 

HEH, 



written as 



fc 
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where h(x, y) is the height of the bilayer relative to a 
flat reference plane, k is the bending rigidity of the 
membrane, and 7 couples the mean curvature to splay 
in m. 
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FIG. 4: Spatial variation of the different order parameters 
in the Pg, phase at T = 310.5 K and C= -4.84 xlO" 19 J. 



FIG. 5: Schematic of the bilayer profile obtained from the 
model. The bilayer thickness is different in the two arms of 
the saw-tooth-like ripples. 



The equilibrium height profile of the bilayer h(x, y) is 
related to the tilt m via the Euler-Lagrange equation, 



V 2 h 



7 



(V-m) 



(4) 



Eliminating h from the free energy density / leads to 
the effective energy density f e ff with a reduced splay 
elastic constant I\ = I\ — j 2 /(2k). 

To determine the mean field phase diagram we choose 
the following ansatz for ip and m, 



ip = ip Q + i/j 1 sin(qx) 



m x = rn . 



m y = mo y 



mix °os (qx) + m 2x sin (qx) 
m 3x cos (2qx) + m^ x sin (2qx) 



(5) 



We do not consider two-dimensionally modulated rip- 
ples, since they do not appear to be generic; there be- 
ing, as far as we know, only one report of such a struc- 
ture (37[. Spatial modulation of m y is neglected as we 



do not keep terms proportional to (V x m) in eqn.© 
for reasons discussed below. Higher order Fourier com- 
ponents of m x are retained in order to account for the 
ripple asymmetry. Three different ripple structures can 
be described with this ansatz: (1) the Pp with no mean 



0), (2) the Pj, with a mean tilt along 
and (3) the P^, with a mean 



tilt (m 0x = m 0y 
V (m 0x = 0, m 0y ^ 0) 
tilt along x (m 0x ^ 0, m 0y = 0). Of these o nly the Pp, 
structure has an asymmetric height profile [33j. 

The phase diagram in the C-T plane obtained from 
numerical minimization of the effective free energy den- 
sity averaged over one spatial period, < / e // > = 

(q/2n) j^^ q feff dx, is given in fig. [3] It is calcu- 
lated for a choice of parameter values, which reproduce 
closely the main- and pre-transition temperatures of 
dipalmitoylphosphatidylcholine at C = - 4.84xlO _19 J. 
The values of T^j a 2 , 03,04, and D are similar to those 



used in refs. [29[ and [34|. As can be seen from fig. 



[3j there are three distinct regions in the phase dia- 
gram corresponding to three different phases: L a (ipo — 
i>i = mox = moy = mi x = m 2x = m 3x = m ix = 0), 
Lp> (ipo 0,m Ox ^ 0, m 0y ^ 0;^i = m lx = m 2x = 
m 3x = m^ = 0), and P/f, (tpo ^ 0,ipi ^ 0,m Ox ^ 
0, m 2x ^ 0, m 3x ^ 0; m lx = m Ax = m 0y = 0). 

It is interesting that of the three possible ripple 
structures only the asymmetric Ph, which is similar 
to the experimentally observed structure, is present 
in this phase diagram [4(|. The first-order transition 
lines which separate these three phases meet at the 
Lifshitz point located at Cl p = —3.32 x 10 _19 J and 
T = 38.75°C. For C > C Lp , the first order L , -> L a 
transition line is parallel to the C axis and occurs at 
T m = 38.75°C. But for C < C Lp , the intermediate P , 
phase is found, separated from the other two phases by 
first-order transition lines. Further, the region occu- 
pied by the Pg< phase expands at the expense of the 
other two as C becomes more negative. 

Typical spatial variation of the order parameters in 
the ripple phase is shown in fig. [4] The height profile 
is asymmetric and resembles very closely those seen in 
experiments (fig. [T]) [l, EH EH . ipi is almost ir/2 out of 
phase with h, so that it is positive (negative) along the 
longer (shorter) arm of the ripple, resulting in different 
bilayer thicknesses in the two arms, again i n ag reement 
with experimental observations (fig. [U 0, lid! ]. Fig. [5] 
shows a schematic of the structure of the bilayer in- 
ferred from these results. It is clear that the model 
presented here accounts for all the salient features of 
the ripple phase listed in the introduction. 

The essential term in the free energy expression re- 
sponsible for the asymmetric ripples is the one with the 
coefficient Tg, since similar structures can be obtained 
by setting r 5 =r 6 =r 7 = 0, as long as r 8 is < 0. Thus 
the present model spontaneously picks out a non-zero 
value of the mean tilt along the rippling direction q, 
even in the absence of any explicit in-plane anisotropy 
of the bending rigidity. This is in contrast to the model 
presented in refs. [20j and [3l|, where the mean tilt oc- 
curs in a direction normal to q, resulting in symmetric 
ripples in the case of achiral bilayers; the bending rigid- 
ity has to be explicitly taken to be lower along the tilt 
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direction in order to obtain a non-zero mean tilt along 
q and to stabilize asymmetric ripples within this model 
[33| . If T-j is made sufficiently positive in the present 
model, ipi and h become almost in phase, so that the 
bilayer thickness is modulated within each arm of the 
ripple. It might be possible to tune this parameter by 
a suitable choice of an impurity which would prefer to 
smoothen variations in ip] in such cases we predict the 
existence of this new ripple phase. 

We have included only terms proportional to (V • m) 
in the expression for / t . In general there will also be 
terms proportional to (V x m), which lead to a rip- 
ple structure with a non-zero winding number in the 



model presented in refs. [20| and [3l|. However, such a 
structure can be expected to be energetically very un- 
favorable in an achiral bilayer, since it is not consistent 
with parallel close-packing of the chains demanded by 
van der Waals interaction. 

A straightforward extension of this model would be 
to use a better description of the chain-melting tran- 
sition, instead of the reduced bilayer thickness ip em- 
ployed here. Possible order parameters include compo- 
nents of the in-plane density wave, as considered in the 
theory of weak crystallization [38j, and those of her- 
ringbone order, used to describe positional ordering in 
monolayers [391 ]. However, such an attempt would be 
useful only if the details of chain ordering in this phase, 
presently unknown, can be experimentally established. 

We thank Yashodhan Hatwalne for many valuable 
discussions and Kheya Sengupta for the electron den- 
sity map of the ripple phase. 
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